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Abstract 

We study the time dependent Hartree equation in the continuum, the semidiscrete, 
and the fully discrete setting. We prove existence- uniqueness, regularity, and approxima- 
tion properties for the respective schemes, and set the stage for a controlled numerical 
computation of delicate nonlinear and nonlocal features of the Hartree dynamics in various 
physical applications. 
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1 Introduction 

In this paper, we study the nonlinear and nonlocal Hartree initial-boundary value problem for 
the (wave) function ■ip{x,t) being defined bj0 

= {-A + v + XV if (x,t) gHx [0,T), 

ij = 0, if (x,t) e 9fi X [0,T), (1) 

ip = ipo, if {x,t) G n X {0}, 

where Q is some domain in M"^ with boundary dQ, and T > is some upper limit of the time 
interval on which we want to study the time evolution of ip. Moreover, v stands for an external 



*aschbacher@ma. tum.de 

^The notation Tp stands for the derivative of ip w.r.t. the time variable t. 
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potential, A denotes the coupling strength, and V is the interaction potential responsible for 
the nonlinear and nonlocal interaction generated by the convolution!^ The system ([T]) has 
many physical applications, in particular for the case Q = W^. As a first application, we 
mention the appearance of within the context of the quantum mechanical description of 
large systems of nonrelativistic bosons in their so-called mean field limit. For the case of a 
local nonlinearity, i.e. V = S, an important application of equation ([T]) lies in the domain of 
Bose-Einstein condensation for repulsive interatomic forces where it governs the condensate 
wave function and is called the Gross-Pitaevskii equation. This dynamical equation, and its 
corresponding energy functional in the stationary case, have been derived rigorously (see for 
example [191 E] ^^^^ [HI; respectively^). In minimizers of this nonlocal Hartree functional 
have been studied in the attractive case, and symmetry breaking has been established for 
sufficiently large coupling. A large coupling phase segregation phenomenon has also been 
rigorously derived for a system of two coupled Hartree equations which are used to describe 
interacting Bose-Einstein condensates (see [HIE] and references therein). Such coupled systems 
also appear in the description of crossing sea states of weakly nonlinear dispersive surface water 
waves in hydrodynamics (see for example [T5l [T8] ). of electromagnetic waves in a Kerr medium 
in nonlinear optics, and in nonlinear plasma physics. Furthermore, we would like to mention 
that equation ([T]) with attractive interaction potential possesses a so-called point particle limit. 
Consider the situation where the initial condition is composed of several interacting Hartree 
minimizers sitting in an external potential which varies slowly on the length scale defined by 
the extension of the minimizers. It turns out that, in a time regime inversely related to this 
scale, the center of mass of each minimizer follows a trajectory which is governed, up to a 
small friction term, by Newton's equation of motion for interacting point particles in the slowly 
varying external potential. Hence, in this limit, the system can be interpreted as the motion of 
interacting extended particles in a shallow external potential and weakly coupled to a dispersive 
environment with which mass and energy can be exchanged through the friction term. This 
allows to describe, and hence to numerically compute, some type of structure formation in 
Newtonian gravity (see [I21II3|)- 

The main content of the present paper consists in setting up the framework for the numerical 
analysis on bounded Q which will be used in [3j for the study through numerical computation 
of such phenomena, like, for example, the dissipation through radiation for a Hartree minimizer 
oscillating in an external confining potential (see also [2j). 

In Section [21 we start off with a brief study of the Hartree initial- value boundary problem ([1]) 
in the continuum setting and we discuss its existence-uniqueness and regularity properties. In 
Section [HI the system ([T]) is discretized in space with the help of Galerkin theory. We derive 
existence-uniqueness and a bound on the L^-approximation error. In the main Section [U we 
proceed to the full discretization of ([T]), more precisely, we discretize the foregoing semidiscrete 
problem in time focusing on two time discretization schemes of Crank-Nicholson type. The first 
is the so-called one-step one-stage Gauss-Legendre Runge-Kutta method which conserves the 
mass of the discretized wave function under the discrete time evolution. The second one is the 



^To be made precise below. 

■^Moreover, the nonlocal Hartree equation has been derived for weakly coupled fermions in [lOj . 
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so-called Delfour-Fortin-Payre scheme which, besides the mass, also conserves the energy of the 
system. We prove existence-uniqueness using contraction methods suitable for implementation 
in [21 [3]. Moreover, we derive a time quadratic accuracy estimate on the L^-error of these 
approximation schemes. In the proofs of these assertions, we write down rather explicit expres- 
sions for the bounds in order to have some qualitative idea how to achieve a good numerical 
control of the fully discrete approximations of the Hartree initial- value boundary problem ([T]) 
for the computation of delicate nonlinear and nonlocal features of the various physical scenarios 
discussed above. 



2 The continuum problem 

As discussed in the Introduction, we start off by briefly studying the Hartree initial-boundary 
value problem ([T]) in a suitable continuum setting. For this purpose, we make the following 
assumptions concerning the domain Q, the external potential v, and the interaction potential V, 
a choice which is motivated by the perspective of the fully discrete problem and the numerical 
analysis dealt with in Section H] and the numerical computations in [3j 

Assumption 1 Q cM.'^ is a bounded domain with smooth boundary dQ. 
Assumption 2 Assumption{J\ holds, v G C^(fi,]R), and V G C^(]R'', M). 
Assumption 3 AssumptionlE holds, and V{—x) = V{x) for all x G M'^. 

The Lebesgue and Sobolev spaces used in the following are always defined over the domain 
fl from Assumption [1] unless something else is stated explicitly. Thus, we suppress fl in the 
notation of these spaces. Moreover, under Assumption [21 let the Hilbert space Ti, the linear 
operator A onTi with domain of definition D{A), and the nonlinear mapping J on be given 



by 



n 

D{A) 



H^nH, 



0; 



A 



-A 



vip + flip] 



where the nonlinear mapping / is defined 



by 



X G 




^Some Hartree-dynamical computations have already been performed in [2]. 



4 



W. H. Aschbacher 



Remark 4 The function A stands for some space depending coupling function which can be 
chosen to be a smooth characteristic function of the domain Q. Such a choice, on one hand, 
insures that all derivatives of f[ip] vanish at the boundary dQ, and, on the other hand, switches 
the nonlocal interaction off in some neighborhood of dQ where, in the numerical computation, 
transparent boundary conditions have to be matched with the outgoing flow of ip (see [2l Ej). 

Remark 5 Under Assumption [21 we have gvW\ ^ C^(]R'^,C) for all ip & L^, and f[ip] G for 
all G L^, see estimate ( ITTj) below. 

We now make the following definition. 

Definition 6 Let Assumption hold, and let T G (0, oo) U {oo}. We call a differentiable 
function ip '■ [0, T) ^ L"^ a continuum solution of the Hartree initial-value problem ([T]) with 
initial condition ipo G D{A) if 

r ii;{t) = Am + vt g [o,t), 

IfT< oo, the solution is called local, and if T = oo, it is called global. 

We make use of the following theorem to prove that there exists a unique global solution of 
the Hartree initial-value problem ([T]) in the sense of Definition [6l In addition, this solution 
has higher regularity properties in time which are required for the bounds on the constants 
appearing in the L^-error estimates in the fully discrete setting of Section HI 

Theorem 7 (Cf. [W, p. 301]) Let Ti be a Hilbert space and A a linear operator on Ti with 
domain of definition D{A) and A* = A. 

(a) Let n E N, and let J be a mapping which satisfies the following conditions for all ip,^ E 
D{A'), 

JD{A^) C D{A^), Vfc = l,...,ri, (3) 

ll^[^]ll < ^^(IIV'II)II^II, (4) 

WA'mW < C{Ul-A\A'-'m\A'n, VA; = l,...,n, (5) 

\\A\.m-jm\ < Vfc = 0,...,rz, (6) 

where each constant C is a monotone increasing and everywhere finite function of all its vari- 
ables. Then, for each ipo G there exists a unique local continuum solution in the sense 
of Definition\^with ipit) G for all t G [0,T). 

(b) Moreover, this solution is global if 

\m)\\<C, Wg[0,T). (7) 

(c) In addition to the conditions in (a), let J satisfy the following conditions for alll < k < n: 
ifijE C^([0,oo),7^) with |r^(t) G for all < I < k and A^^-^^^jit) G C{[0,oo),n) 
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for alio < I < k, then 

J[?/'(t)] is k times dijferentiable, (8) 
^Jim] e D{A--'^-'), (9) 

J k 

A--'-'—J[m] e Ci[0,oo),n). (10) 

If this condition holds, then the local solution from (a) is n times differentiahle and -^il){t) G 
/^(A"-^) for alll<k<n. 

Remark 8 With the help of Theorem [7] (and its proof) the constants in the estimates fl7^ and 
( I87p below on the L^-error of the fully discrete approximations are finite and can be estimated 
explicitly. 

Proof Let us start off by checking condition (jlj). 

Condition ©: H^ML^ < C(||V^||^.)II^IIl^ 
Using Assumption [2] and the estimate 



l|A^y[<^^]x|lL2 < ||A||^oc||1^1|ioo(Kd)||^|lL2||^/'l|i2l|xllL2, 

we immediately get 



'111 



and the pref actor C(||'?/'||) from (jl]) is a monotone increasing function of ||'?/'||^2- Next, let us 
check condition 

Condition \\A^J[i)]\\^2 < C(||^||^2) || A^ILa for all A; e N 

In order to show we have to control the L^-norm of A''{vip) and of A^(A(7y [iV^p]-?/^). Hence, 
we write the powers of the Dirichlet-Laplacian as follows, 

«i,/3i,...,Ofe,/3fe=0,l 

^\9v[m =E c^.P...aA9,^.,,, ,^,g0,y[v]id^r'9]^^^^^^ (13) 

, Jl Jl'" Jfc Jfe 

JivJfe=lv,a 
oi,/3i,...,Ofe,/3fe=0,l 



where Cai/3i...aki3k denote some combinatorial constants. Hence, using ffTTj) and the following 
Schauder type estimatelf] 



'I.e. the elliptic regularity estimate of generalized solutions up to the boundary, see for example [521 p.383]. 
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we get the following bounds on ffT^ and ffT^ . 

w^'i^niL^ < CM 

< C||^||^^...^||A^^||^„ (14) 
Using ( |T^ . and ( |T5i) twice, we arrive at 

)I|a'=^1Il2, (16) 

where the prefactor CdlV'H, H^'^ll, •••5 1|^'^~^'^||) from depends on ||'?/'||2,2 only and is mono- 
tone increasing in ||'?/'||^2. Next, we check condition ([3]). 

Condition (E]): JD{A^) C D{A'') for all G N 

Due to (UnD and since D{A'') = {ijj e H'^ n \ AiP e n Hi, A^^^V e if^ p iJi}, it 
remains to show that A^ J{iIj) G for all ip G D{A'') and for all j = 0, k. But this follows 
since f , A G and from the fact that C°°{Q) is dense in if™ w.r.t. the ii^-norm for all 
m G Nq. Let us next check condition ([U]). 

Condition ®: ||A'=(J[^] - J[e])||^2 < ^^(11^1^2, UWl^, IIA^IIl^, || A^^ILO || A^ - A'=eiL2 
In order to show we write the difference with the help of the decomposition 

+ fi'yb2^2](xi - X2)- (17) 

Each term on the r.h.s. of ffTTj) can then be estimated with the help of ffT^ . Hence, as in 
we get 

||A'=(J[^]- J[e])||^2 < Ci\\v\\^,,^^ + \\\\\ 

■(ll^lli^ + (11^11,2 + ||e||,2)||AVlL2)) ||AV - A^^ll,., (18) 

where the prefactor Cdl'i/'ll, ||^||, ||A'^'?/^||, HA'^^H) from ([6j) depends on the lowest and the 
highest power of A only, and it is monotone increasing in all its variables. 

Condition dH): ||^(t)||i2 = H^oIIl^ for all t G [0,T) 
This condition is satisfied due to Proposition [H] below. 

Condition ([H]): J[ip{t)] is k times differentiable in w.r.t. time t 

Let n G N be fixed, and let k = 1. Then, it is shown in [IK, p. 299] that part (a) and 
conditions (IHD, and ([10]) for A; = 1 imply that ^ G C^{[0, 00), L^) with ^V(^) e D{A''~^) 
and v4"-2^V(i) e C([0,oo),L2). Then, using conditions ([H]), ([H]), and ffTUl) for subsequent 
1 < < n leads to the claim of part (c) by iteration. Hence, we have to verify that conditions 
(IH]), Q, and ( ITOi) are satisfied for 1 < k < n. To this end, we make use of decomposition ( |T7I) 
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to exemplify the case k = 1 and to note that the cases for /c > 2 are analogous. In order to 
show that J[ip{t)] is differentiable in L^, we write, using ffTTl) . 

J[^(t + h)] - J[^(t)] ^ ^{t + h)- ^(t) 
h ~ ^ /i 

+ Xgy[ij(t + h)i^i±^^]ij{t + h) 

Applying ( ITTl) and ( fTTl) . we find that J[ilj{t)] is differentiable in w.r.t. t with derivative 

^j[^(t)] = + A^7y[^(t)^(t)]^(t) + \gv[mm]m + \gv[\m\']m- m 

Making use of (dH), (dH), and (USD in the estimate of || A"-2(^ _ ^ similarly 

to (HHD, we find that G D(A"-2) and that A'^-^^ J[^(t)] G C([0, oo), L^). Due to the 

structure of ( fT9l) . we can iterate the foregoing procedure to arrive at the assertion. □ 



In order to verify condition ([7]), we define the mass A^f?/'] and energy T-C[ip] of a function G 
by 

M[^] := ml., 

We then have the following. 

Proposition 9 Let Assumption hold, and let ip he the unique local continuum solution of 
Theorem^ Then, the mass and the energy of ip are conserved under the time evolution, 

M[iP{t)] = M[iPol VtG[0,T), (20) 
nim] = ^[^o], VtG[0,T). 



Proof From Theorem [7]it follows that the function t Ai[ilj{t)] belongs to C^{[0,T),Wq) with 

^A<[^(t)] = 2Re(^(t),^(t))^„ 

which vanishes due to ([2]). For the conservation of the energy, we have the following three 
parts. First, using the regularity of ip in time and —Aip)^. = WVipWl. for all ip G H"^ fl Hq, 
we observe that t ^— > ||VV^(t)||^2 belongs to C"^([0, T), Rq ) and has the derivative 

-^||V^(t)||i. =2Re(^(t),-A^(t))^.. (21) 
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Second, the function t i— > {ip{t),vip{t))^2 belongs to C"'^([0, T), M) and has the derivative 

^^{ip{t),vm)L^ = '^^eim,vm)L^- (22) 

Third, using I'j/'p— = ip{ilj—i^) + {ip—ip)ip in the decomposition of /[■?/']—/[(/)] as in (fTTIl . we get 
^^y = gv[\'ip\'^]^ + 2Re igv[i>ij]) ^ in L^, and therefore the function t (ipit), /[^/'(t)])i2 

belongs to C^([0,T),M) and has derivative 

^(^(t),/[^(t)])^, =4Re(^(t),/[z^(t)])^„ (23) 

where we used Assumption [3] to write {g[ijjtp]'ilj , ■?/')^2 = {ip, /['^])l2. Finally, if we take the scalar 
product of (I2l) with 'ip and the real part of the resulting equation, we get 

= Re iil^WIll^ = Re [(^(t), -A^(t))^2 + {m,vm)L^ + (^(t), /[^(t)])^^] • (24) 
Plugging (I2T1) . (!22l) . and (!23l) into (Eij), we find the conservation of the energy H['ilj{t)]. □ 



Remark 10 For more general interaction potentials V, in particular in the local case f[%l)\ = 
l^l'^ip in = 20 one can use an estimate from [7j which controls the L°°-norm of a function 
ip G by the square root of the logarithmic growth of the if^-norm. 



<C(l + y/log(l + 11^11^2)). 



where the constant C depends on HV'lliyi- This estimate allows to bound the graph norm of the 
continuum solution by a double exponential growth, and, hence, makes the solution global. 

Taking the L^-scalar product of w.r.t. functions ip G H^, and using again that {ip, —Aip)^2 = 
(Vif, Vip)i^2 for all ip G fl Hq and for all ip G Hq, we get the following weak formulation of 
the continuum problem ([21), 

i(¥P,^)^2 = (Vyp, VV^)^2 + i^,viP)^2 + i^,m)L2, Vy. G H^, Vt G [0,T), , . 

m = v^o. ^ 

The formulation ( l25i) is the starting point for a suitable discretization in space of the original 
continuum problem. We will discuss such a semidiscrete approximation in Section [3l 



3 The semidiscrete approximation 

In this section, we discretize the problem fl^^ in space with the help of Galerkin theory which 
makes use of a family {Sh}h£(o,i) of finite dimensional subspaces approximating the infinite 
dimensional problem in the following precise sense. 



^We are mainly interested in d = 2 for the numerical computations. 
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Assumption 11 The family {Sh}he{o,i) of sub spaces of Hq has the property 

Sh C cm n H^, dim St = Nh < oo, V/i G (0, 1). 

Remark 12 For tlie numerical computation in [3], the physical space is (a smoothly bounded 
superset of) the open square VL = (0, D)"^ C with D > whose closure is the union of the 
{n — 1)^ congruent closed subsquares generated by dividing each side of fl equidistantly into 
n — 1 intervals. Let us denote by = {n — 2)^ the total number of interior vertices of this 
lattice and hy h = D/{n — 1) the lattice spacing0 Moreover, let us choose the Galerkin space 
Sh to be spanned by the bilinear Lagrange rectangle finite elements (pj G C{Q) whose reference 
basis function ipo : Q ^ [0, oo) is defined on its support [0, 2h]^^ by 



xy, if {x, y) G [0, h] 



x2 



^^,A._1J {'^h-x)y, if(x,y)G[/i,2/.]x[0,/i], 

V^{x,y).-^< ^2h-x){2h-y), if {x,y) e [h,2h]^\ ^^^^ 
[ xi2h-y), if{x,y) G [0, /i] x [h,2h], 

see Figure [H The functions (fj are then defined to be of the form ( l26i) having their support 
translated by (rnih,m2h) with mi,m2 = 0, ...,m — 3. Hence, with this choice, we have C 
cm n H^{n) and dim Sj, = N^. 




Figure 1: ipo{x,y) on its support [0,2/i]^^ with maximum at vertex {h,h). 



Motivated by the weak formulation (!25l) . we make the following definition. 

Definition 13 Let Assumptions andlIJ\ hold. We call iph '■ [0,T) Sh with iphy'^Ph ^ 
L^(0,T; Sh) a semidiscrete solution of the Hartree initial boundary-value problem ([T]) with initial 
condition ipoh ^ Sh if 

r iA(y,,^^)^^ = {Vip,ViJh)L^ + {^,v^h)L^ + {ipjm)L2, e Sh, Vt G [0,T), 

^As bijection from the one-dimensional to the two-dimensional lattice numbering, we may use the mapping 
r : {0, ...,TO — 1}^^ — > {0, ...,m^ — 1} with j = r(mi,m2) := mi -I- m2m. 
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Remark I4 In general, the weak problem (1231) is set up using the Gelfand evolution triple Hq C 
L2 c (H^Y = H-\ One then looks for weak solutions iIj e W^iQ.T] Hi, L"^) C C([0,T),L2) 
motivating Definition [131 

We assume the Galerkin subspace Sh from Assumption [TT] to satisfy the following additional 
approximation and inverse inequalities. 

Assumption 15 Let Assumption [711 hold. Then, there exists a constant Ca > s.t. 
inf (11^ - ^\\^, + hU - vWm) < CAh'Mn^, ^ijeH'n Hi 

Remark 16 For an order of accuracy r > 2 of the family {Sh}h&{o,i)-, the usual assumption 
replaces the r.h.s. by C^/i'^ll'^llji^s and is asked to hold for all ip G H^ fl H^. For simplicity, we 
stick to Assumption [T5l 

Assumption 17 Let Assumption [771 hold. Then, there exists a constant Cb > s.t. 

Ilv^ll/^i < C's/^^'^ll'y^ll^a, e Sh- 

Remark 18 For the two-dimensional bilinear Lagrange finite element setting of Remark [121 
both Assumption [TJl and Assumption [T7| holdj^ 

Furthermore, we make an assumption on the approximation quality of the initial condition 
V'oh G Sh of the semidiscrete problem ( 1271) compared to the initial condition ipQ e H'^ fl H^ of 
the continuum problem Q- 

Assumption 19 Let Assumption [771 hold. Then, there exists a constant Cq > s.t. 

11^0-^0/^11^2 < Co/i^ (28) 

The semidiscrete scheme has the following conservation properties. 

Proposition 20 Let Assumptions l3l and 1771 hold, and let iph be a semidiscrete solution of the 
Hartree initial boundary-value problem (til) in the sense of Definition [73 Then, the mass and 
energy of iph are conserved under the time evolution, 

M[Mt)] = MiijQhl Vt G [0,r), (29) 

Proof If we plug ip = iphit) into ( 1271) and take the imaginary part of the resulting equation, we 
get the conservation of the mass. If we plug p = iphit) into ( 1271) and take the real part of the 
resulting equation, we get the conservation of the energy using Assumption ISlI^ □ 

Existence-uniqueness is addressed in the following. 

^See for example P", p. 109, 111]. 
^In the sense of Remark HH 
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Theorem 21 Let Assumptions [3l and [77] hold. Then, there exists a unique global semidiscrete 
solution iph of the Hartree initial- boundary value problem ([T]) in the sense of Definition[ 



Proof Let {ipj}^J^^ be a basis of the Galerkin space Sh, and let us write 

i=i 

Plugging (IHUj) into the semidiscrete system flTTj) . we get for z(t) := {zi(t), zjSf^it)) G C^**, 

izit) = A-\B + Y)z{t) + A-^H[z{t)]z{t), Vt G [0, T), 
z(0) = 2:0, 

where V'o/i = Tl^=i{^o)j '■Pj ^^nd the matrices A, B E C^^xn^ g^j.^ ^-j^g positive definite mass and 
stiffness matrices, respectively, 

Aij := {(pi,(pj)^2, 
Bij := (yipi,V(pj)^2- 

Moreover, Y G C^'"^^'' is the external potential matrix, 

and the matrix- valued function H : C^'" QNhxNh jg (defined by 

k,l=l 

Since the function C^*" 3 z ^ A~^{B + Y)z + A~^H[z]z G C^*" is locally Lipschitz continuous 
analogously to the continuum case, the Picard-Lindelof theory for ordinary differential equations 
implies local existence and uniqueness of the initial value problem ( l3T|) . Moreover, this local 
solution is a global solution if it remains restricted to a compact subset of C^'" . But this is the 
case due to the mass conservation from fl29l). □ 



We next turn to the L^-error estimate of the semidiscretization. For that purpose, we introduce 
the Ritz projection^ 

Definition 22 Let Assumption [77] hold. Then, the Ritz projection : Hq Sh is defined to 
be the orthogonal projection from onto Sh w.r.t. the Dirichlet scalar product (V-, V-)^2 on 
Hi, i.e. 

(V(/^,Vi?,,^)^2 = (V¥^,V^)^„ ^ipeSh. (32) 

^°Also called elliptic projection. 
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The Ritz projection satisfies the following error estimate. 

Lemma 23 (Cf. 120, p-8]) Let Assumptions [77] and [75l hold. Then, there exists a constant 
Cn>0 s.t. 

11(1 - RM\^, + hUl - RnmH^ < Cnh'UW^,, e n Hi (33) 
The next theorem is the main assertion of this section. 

Theorem 24 Let Assumptions [77| [75|, and\T^hold, and let ip be the solution of the contin- 
uum problem from Theorem and iph the solution of the semidiscrete problem from Theorem 
[m Then, for any < T < oo, there exists a constant Ce > s.t. 

max llV'(t) -^^/^(^)llL2 < CeH^ 
Proof We decompose the difference of and iph as 

ij{t)-Mt) = pit)+m, 

where p{t) and 6{t) are defined with the help of the Ritz projection Rh from (!32i) by 

p{t) := (1 - Rh)m, d{t) ■= RhiPit) - iPh{t). 
Making use of the schemes fl25l) . fl27|) . and (|32|) . we can write 

i(<^, e{t))^, - (v^, v^^(t))^. = -i(<^, p(t))^. + (<^, v{ip{t) - ^,(t)))^. 

+ (v9,/[^(t)]-/[V^,(t)])^.. (34) 

Plugging if = 0{t) E Sh into ( l34l) and taking the imaginary part of the resulting equation, we 
get the differential inequality 

\ ^ \m\\% < (iipwiIl^ + \\vm) - i^Hmy + wmm - fiumy) 

< mt)\y+c^mt)\y + \\e{t)\y))\\^{t)\y. (35) 

where we used the conservation laws (^01) and (^Uj) and (|TS|) to define the constant 

•= II'^IIl°° + 2 II a II ^oo II ^11 Loo(]f5d) 

;ularize the 

asi;|ll^lli. = |^(ll^lli.+a,weget 



Using e > to regularize the time derivative of ||^||^2 at = by rewriting the l.h.s. of ([35 

1 d_||/]||2 _ l_d 



^ (ll^lli. + e')'" < Wmh^ + (llp(t)IL. + ||^(t)IL.) , (36) 
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1 /2 

where we used ||j;^2 < (II^ILa + e^) • Integrating (|36i) from to t, letting e — * 0, and 
applying Gronwall's lemma to the resulting inequality, we find 



mt)\\L^ < ii^(o)iL.+ / d.(iip(.)ii,.+ciiip(.)ii,o 

Jo 

+ c,j\s (\\emy + [^n (IIpHII,. + c,\\p{u)\y)) e-(*-^). (37) 

In order to extract the factor /i^, we apply (|33l) and Assumption [T9] to get 

||p(t)IL. < Cn\\m\\H^h\ (38) 
||p(t)IL. < Cn\m)\\H^h\ (39) 
11^(0)11^. < {C, + CR\m^,)h\ (40) 

Plugging ([3H1), dMD, and (HOD into dSZD, we finally arrive at 

wm-umL^ < Pit)\\^. + iip(t)ii^. < c,{t)h\ 

where the time dependent prefactor is defined by 



C2{t) := {Co + Cn\\MH^y'' + CR ds{ms)\\H,+cMs)y,) 

Jo 

+c,Cr [ ds rd«(||^(n)||^.+ci||^(w)||^.)e'^^*"'^ + ^«II^WII//- 



^0 JO 

Setting Ce '■= niaxtg[o,T] C2(t) brings the proof of Theorem [2^ to an end. □ 



Remark 25 For the local case f\%l)\ = one replaces the original locally Lipschitz non- 

linearity / by a globally Lipschitz continuous nonlinearity which coincides with / in a given 
neighborhood of the solution ip of the continuum problem. One then first shows that the 
semidiscrete solution of the modified problem satisfies the desired L^-error bound, and, second, 
that for h sufficiently small, the modified solution lies in the given neighborhood of ip. But for 
such h, the solution of the modified problem coincides with the solution of the original problem, 
and, hence, the solution of the original problem satisfies the desired L^-error bound, too. 



4 The fully discrete approximation 

In this section, we discretize the semidiscrete problem ( 1271) in time. To this end, let us denote 
by e N the desired fineness of the time discretization with time discretization scale r and 
its multiples tn for all n = 0, 1, 2, A^, 



T 

T := — , tn := nr. 



(41) 
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As mentioned in the Introduction, we will use two different time discretization schemes of 
Crank-Nicholson type to approximate the semidiscrete solution i/j^ of Theorem [21] at time t„ 
by \E'„ G \E', where 

^■.= {^o,^i,...,^N)es;:^''^'l (42) 

These two schemes differ in the way of approximating the nonlinear term follows. 
Let Af := {1, 2, N} and Afo := Af U {0}, and define 

^„-i/2 := ^ (^n + , G AT. (43) 

Then, the first scheme implements the one-step one-stage Gauss-Legendre Runge-Kutta method 
in which the nonlinear term is discretized by 

9v[\^n-l/2\^]-^n-l/2- (44) 

In this method, the mass A^[\I'„] is conserved under the discrete time evolution. The second 
scheme, introduced in [S] and applied in [T], discretizes the nonlinear term by 

(7v[|(|^„|' + |^n-in]^n-l/2. (45) 

This method, in addition to the mass, also conserves the energy 7i[\E'„] of the system. In the 
following, for convenience, we will call the first scheme coherent and the second one incoherent. 



Coherent scheme 

In order to define what we mean by a coherent solution of the Hartree initial-boundary value 
problem ([T]), we define 

■■= - (^n - ^n-l) , Vn G A^. 
T 

Definition 26 Let Assumption IE and fill hold. We call ^ G S^^^^^'' a coherent fully discrete 
solution of the Hartree initial-boundary value problem ([T]) with initial condition ipQ^ G Sh if 

G Sh, Vn G AT, (46) 

^0 = i^oh- 

The coherent solution has the following conservation property. 

Proposition 27 Let \E' G S^^^'^^^ be a coherent fully discrete solution of the Hartree initial- 
boundary value problem ([1]) in the sense of Definition IME Then, the mass of \I' is conserved 
under the discrete time evolution, 

M[^n]=M[iJoh], VnG A/'o. (47) 
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Proof If we plug ip = into fHUl) and take the imaginary part of the resulting equation, 

we get 

= Imi(^„_i/2, = ^ {M[^n] - A^[^„-i]) . 

□ 

Remark 28 The energy 7i[\l/„] of the coherent solution fl46l) is not conserved under the discrete 
time evolution, see [T] and references therein, in particular [17] and |21J for the local case with 
d=l. 

The question of existence and uniqueness of a coherent solution is addressed in the following. 
Theorem 29 Let Assumptions [H [77], and [7^ hold, and let the time discretization scale r he 



sufficiently small. Then, there exists a unique coherent fully discrete solution of the Hartree 
initial-boundary value problem ([T]) in the sense of Definition\ 



Proof Let G be given, and define the mapping : Sh Sh by 

{^,F4ij])^, := - ^ ((Vy., VV^)^, + + iy^jmL^), Vy. e S,. (48) 

For some n G A/", let the n-th component ^^n-i of \E' G S^^^'^^^ from fl42p be given. Adding 
2i [ip, \E'„_i)/r on both sides of ( l46ll . we can rewrite (l46l) with the help of (HSl) in the form of a 
fixed point equation for 

vI/„_i/2 = F*_Jv&„_i/2], (49) 

from which we retrieve the unknown component \Ef„ by 0431) . In order to construct the unique 
solution of ( H9I) . we make use of Banach's fixed point theorem on the compact ball Bn-i '■ = 
{ip e Sh \ \\ij\\L2 < A^[^„-i]^/^ + l} in Sh. Using Assumption [IT] and 1^, we get, for ip,^ e Sh, 



\{^,F^„_,[ij]-F^„_M 



L2\ 



Plugging 99 = F^ij^_-^^[ip] — into flHJ|) and picking ip and ^ from Bn-i, we find 

\\F<,^-A^]-F^^_AmL^ < (51) 

where the constant q;„_i is defined, for all n G Af, by 

:= Clh-^ + \\v\\^^ + 4||A||^^ II^IIloo(k.)(-M[^„-i]'/' + 1)^ (52) 



Let now an-i{M[ifn-iY^^ + l)r < 1. Then, it follows from (EI]) and ([52]) that Fvi>„_j maps Bn-i 
into Bn-i (set ^ = in ( l5l]) ) and that -Fii-^.i is a strict contraction on Therefore, for 
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such r, Banach's fixed point theorem imphes the existence of a unique solution '^n-1/2 £ -Bn-i 
of the fixed point equation fH^ . Moreover, due to the mass conservation (jUj), there exists no 
solution \E'„_i/2 of ( l49l) with \E'ri--i/2 & Sh \ Bn^i- Hence, the component of the coherent 
solution exists and is unique for such r. Starting at \E'o = i'oh and proceeding iteratively, we 
get all n + 1 components of the coherent solution G S^^^~^^\ Moreover, again due to ( H71) . 
we get a uniform bound on the size of the time discretization scale r, e.g. 

ao(-M[^oh]'/' + l)r<l. 

□ 

Remark 30 Since ao > C'^h^'^, we have that r < C^^/i^, where Cb stems from Assumption [T71 

We next turn to the first of the two main assertions of the present paper which is the time 
quadratic accuracy estimate on the L^-error of the coherent solution. 

Theorem 31 Let Assumptions [771 [75l andUM hold, and let G S^^'^'^^^ be the coherent 
solution from Theorem\2^ Then, there exists a constant Ck > s.t. 

max \\^{tn) - "^uWl^ < Ck{t^ + h^). 

n£Mo 

Remark 32 The constant Ck depends on higher Sobolev norms of the continuum solution ip. 
These norms exist due to the regularity assertion in Theorem [7] (c). 

Proof Let G A/" be fixed and define ipn '■= V'(^n) with t„ from fj^Tl) . As in the proof of Theorem 
[Ml we decompose the difference to be estimated as 



V'n - ^n. = Pn + (53) 

where p„ and On are again defined with the help of the Ritz projection from ( [321) by 

Pn := (1 - Rh)'(pn, On := Rh^n " ^n- (54) 

Using Taylor's theorem in order to expand ipn around t = up to zeroth order in tn and the 
estimate on the Ritz projection ( l33i) . we immediately get 

\\Pnh2<CRh^{\mH^ + l^^dt ||^(t)||^.). (55) 
In order to estimate 0^, we want to extract suitable small differences from the expression 

Ln.,fi ■= -{'^,0n - On-l)L2 - ^(Vv?, V(6'„ + 6'„,-i))^2 - ]-{ip,v{On + On-l)) ^2 (56) 

r z z 



which contains all the linear terms in (1461) moved to the l.h.s. with \E'„ replaced by On- For 
this purpose, we first plug the definition of On into f[3UI) . and then use the definition of the Ritz 
projection f[H^ and the scheme (jlUI) to get 



Ln,^ = ^{ip,Rh{i^n-i^n-l))L2 " ^ ( Vv^, V(^n + ^n-l ) )l2 - ^{ip, vRhH^n + i^n-l)) 

-(¥.,/[vl>„_i/2])i2. (57) 
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Rewriting the first term on the r.h.s. of fl37|) with the help of the continuum solution satisfying 
the weak formulation (1^. we have 



L2 



- (V3, {Rh - 1)(^„ - V'n-l))L2 + i (V5, ^{Ai - ^n-l) " ^n-l/2)i2 
+ (V<^, V^„_i/2)^2 + (<^, ^^V'n-l/2)i2 + (<^, /[^n-l/2])i2, (58) 



where we used the notations 4'n-i/2 '■= i'itn — t/2) and '■= i^(tn — ''"/2)- Plugging 

into (1H7|) . we can express in the form 



L2' 



(59) 



f 1) 

where the functions oon with j = 1, 6 are defined by 



(1) 

n 

(2) 
n 

(3) 
n 

(4) 
n 

(5) 
« 

(6) 



-{Rh - l)(V'n - V'n-l), 

r 

i(7(V^n - V^n-l) -^n-1/2), 
A(i {A + tpn-l) - V'n-1/2) , 
V {lpn-l/2 - K^n + V'n-l)) , 

Rh) {ipn + V'n-l) , 



fbPn-1/2] - /[*n-l/2], 

and, for ujn \ we used again (Vy^, V'?/')^2 = (v^, — A'?/')^2 for all G i?^ fl Hq. Plugging ip = 
{On + ^n-i)/2 into and (13^ . and taking the imaginary part of the resulting equation, we 
gefl 

6 



l^nll 2,2 < 11^2 + T 



(i)| 



Il2- 



(60) 



i=i 



Let us next estimate the terms ||a;n ■'||j;^2 for all j = 1,...,6. For u;n \ we expand ipn around 
t = tn-i up to zeroth order in r and use ( l33i) s.t. 



(1)1 



^2 ^ L'iJ 



2_-l 



dt ||^(t)||^2. 



(61) 



tn-l 



(2) 

For Un , we expand V^n-i and ipn around t = tn — t/2 up to second order in r/2. 



' " - 2r 



< 



tn-T/2 



dt {tn-l-tyUm^, + / dt (t„-t) 



tn-T/2 



L2 



dt 



L2- 



(62) 



^If On = 0, we are left with jSSl). 



18 



W. H. Aschbacher 



Analogously, for Un^ and ujn \ we expand ipn-i and ipn around t = t„ — r/2 up to first order in 
r/2, 

< I WArnh^, (63) 

< ^II^IIlo^ / dt (64) 

^ in — 1 

(5) 

For ujn , expanding V'n-i and ifjn around t = up to zeroth order in time, we get, analogously 
,(1) 



to the estimate of ufi , 



\W^h^ < Cnh'\\v\\^^(^\\Mm + I \\m\\H^)- (65) 
Finally, for ujn \ we apply the local Lipschitz continuity ffTSl) to get 

||t^i^^|lL2 < Ci ||V'„_l/2 - ^n-l/2|li2, (66) 

where we used the continuum mass conservation (1201) and the coherent fully discrete mass 
conservation fHTl) to define the constant 



ci ■.= 2\\\\\^^\\V\\^^^^,^{M[iJo]+M[iJoh]). (67) 
Since we want to reinsert the decomposition ( l53l) into the r.h.s. of ( l66l) . we write 



IIV'n-l/2 - ^n-l/2||j^2 < ||^n-l/2 " |(V^n + V^n- 



+ ^(l|Pn-l|lL2 + IIPnlL^ + ll^n-l|lL2 + ||^n||i2). (68) 



Plugging the estimates (13^ . (lUTj) to (lUUj) . and (lUHj) into (p7|) . we find 

Cl 

2 



enh.< WOn-ih. + (A« + TA^^^y + + ^ r (||^^„_J,. + ||^^„||^.) , (69) 



where the first term on the r.h.s. of flU^ was estimated as in uJn^ or ujii'\ and 

:= Cn / dt ||^(t)||^2, (70) 
:= C«(ci + ||t;||^^)(||^o||^,2+ Tdt ||^(t)||^.), (71) 





1 rtn 1 /"in 1 rtn 

:= 7(ci + II^IIlo=) / dt||^(t)|L. + - / dt||^'(t)||^.+ - / dt ||A^(t)||^.. (72) 



^ ^n — \ ^ i^n — 1 ^ t-Ti — 1 

If we choose the time discretization scale r to be small enough, e.g. C\T < 1, we can construct 
the following recursive bound on ||6'„||2,2 from inequality flU^ . 

||^n|L2 <i?(^)||^^„-l|L2+Bf , (73) 
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where we used that 1/(1 — Cit/2) < 1 + Cit if cit < 1 to define 

B^'^ := (1 + Cirf, 

5f := (l + Cir)((A«+r4^))/.^ + A[f)r^). 
Therefore, if we iterate the bound (ff3|) until we arrive at ||6'o||^2, we get 

n 
k=l 

n 

< c.(||^olL. + E((4^^+-4'^)^' + 4V)), (74) 

k=l 

where, on the second fine of ( 1741) . we first extract the global factor (B^^^)" which can then be 
estimated as (B^^))" < (1 + CiT/Nf^ < C2 with the definition 

(75) 

It remains to estimate ||^o|Il2 on the r.h.s. of (17^ . This is again done by using the estimate 
on the Ritz projection (1331) . 

\\0o\\l2 < 11^0 - V^o/.|Il2 + II (/?/,- 1)^01^2 

< Uo-i^oh\\L2+CRh^\\MH2- (76) 

Hence, with estimate (13^ on ||pn|l2,2 and the estimates (!7il) and (!7U|) on ||6'„||^2 in the decom- 
position (!54l) . taking the maximum over all times, we finally arrive at 

max \\lpn - "^nWii < C2\\i'0 - V'0/i|Il2 + Cs/l^ + C4T^, (77) 

where the constants C3 and C4 are defined by 

C3 := CK(l + C2(l + T(ci + ||t;||^^)))(||V'o|lH^ + |'^dt ||^(t)||^2), (78) 

T 1 pT 



C4 := I ((ci + ||.;|L.) 1^ dt Wmh^ + 2 io ""^'^^^"^^ ^ /o "^"^^^^"^O- (^^^ 



The constants Ci and C2 are given in flU7|) and fl73|) . respectively. Using Assumption [T^ and 
setting := max{c2Co + C3, C4} brings the proof of Theorem [31] to an end. □ 



Remark 33 We can compute an explicit bound on the integrands in ( ITHll and (1791) on any finite 
time interval. As an example, for the last term in the constant C4 from (1791) . we havj^ 



i^See Theorem [71(c) and HSl p.299]. 
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The first term is exponentially bounded in time using the conditions from Theorem [7] (a) and 
(b) and Gronwall's lemma on the Duhamel integral form of the differential equation ([2]), 

where the constant C stems from f|T6l)l^ The second term is again bounded due to ([5]). Finally, 
the third term is bounded due to equation ffT^ for the time derivative of the nonlinear term 
J[ip{t)] and the corresponding estimates ffT^ and ffT^ . 

Incoherent scheme 

As described at the beginning of Section HJ we also study a second discretization scheme which 
approximizes the nonlinear term (yfy not by (jSj) but rather by the expression ( H5l) . 

Definition 34 Let Assumptions IE andlll\ hold. We call \1' G S'^^'^^^^ an incoherent fully 
discrete solution of the Hartree initial-boundary value problem with initial condition ipQh G 
Sh if 

Vv? G Sh, G Af, (80) 

^0 = i'oh- 

The incoherent solution has the following conservation properties. 

Proposition 35 Let \E' G S*^^^^^"* be an incoherent fully discrete solution of the Hartree initial- 
boundary value problem in the sense of Definition and let Assumption hold. Then, 
the mass and the energy of \E' are conserved under the discrete time evolution, 

M[^n] = M[^oh], VnGA^o, (81) 

Proof Plugging ip = '^n~i/2 into ( IHOl) and taking the imaginary part of the resulting equation 
leads to the mass conservation as in the proof of Proposition [271 In order to prove the energy 
conservation, we plug = \E'„ into (IHOj) and take the real part of the resulting equation. Using 
that (^'„,A^y[|^^„_i|2]^'„)^2 = (^'„_i,A^y[|^„|2]^„_i)^2 duc to Assumption El we get 

= Re i(^„, ^„)^2 = ^ m^n] - n[^n-i]) . 

IT 

□ 

We next turn to the proof of existence-uniqueness of the incoherent solution. 



^•^See also [TO, p. 300]. In contradistinction to the general case from Theorem[7] (a), the growth rate C from 
pop only depends on the mass of the initial condition (and on w, V ^ and A, of course). 
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Theorem 36 Let Assumptions [77], and [7^ hold, and let the time discretization scale t be 
sufficiently small. Then, there exists a unique incoherent fully discrete solution of the Hartree 
initial- boundary value problem ([T]) in the sense of Definition l34\ 

Proof The proof for the incoherent solution is analogous to the proof of the coherent solution. 
Let 4) E Shhe given, and define the mapping : Sh ^ Sh by 



L2 



~j ((V^, V^)^. + {if, v^)^, + (yp, \Xgv[\2^lj - 0p + 101^)^0 , G Sk. (82) 

For some n G A/", let the n-th component ^n~i of G S^^^^^^ from ( H2l) be given. Adding 
2i((y9, \E'„_i)/r on both sides of flHOl) . we rewrite (IHOj) with the help of fl82|) in the form of a fixed 
point equation for \E'„_i/2, 

In order to make use of Banach's fixed point theorem as in the proof of Theorem [221 we show 
that maps the compact Ball Bn-i '■= {if) E Sh \ HV'IIlz < A^[\E'n-i]"'^''^ + 1} into itself and 

that is a strict contraction on Bn-i- To this end, we write 

<^(l|V(^-e)llL2||Vv9||^. + l|t;||^^||^-eiL.||¥^L.+ ^AIIV'IL.), (83) 

where, with the help of ( fTTl) . ( fT71) . and \ \z\'^ — < + — w| for all G C, the third 
term A on the r.h.s. of (l83l) can be estimated as 



A := ||A^v'[|2^-^„-iP + |^n-i|']^-A(7v[|2e-vI>„-iP + |^„^iP]eL2 



Hence, plugging ip = G*„_J'?/'] — into ( !83|) . we get for ip,^ E Bn-i using Assumption 

where := + HfH^cx. + 12||A||^oo||^||ioo(iRd) (A/li^E'n-i]^/^ + 1)^ like in the coherent 

scheme (l52l) . Therefore, we arrive at the claim as in the proof of Proposition [29] using the mass 
conservation from flHTl) . i.e. the incoherent solution exists and is unique if the time discretization 
scale T is sufficiently small, e.g. ao{A4[ipohY^'^ + l)r < 1. □ 



Finally, we also provide a time quadratic accuracy estimate on the L^-error of the incoherent 
solution. Again, the proof is analogous to the corresponding proof for the coherent solution 
from Theorem [3TJ 
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Theorem 37 Let Assumptions QT], [731 andUM hold, and let \E' G S^^^~^^^ be the incoherent 
solution from TheoremlEE Then, there exists a constant Cj > s.t. 

max\\i;{tn) - ^n\\L2 < Ciir^ + h^). 

Proof As in f l3H|) and fl3^ . we make use of the decomposition ^/^^ — \E'„ = pn + On, and we 
estimate again by ( l55l) . In order to estimate ||6'n||j;^2, we define Ln,ip as in ( l56l) . Then, 
everything in equation ( 1571) remains unchanged up to the last term which is replaced by the 
expression —{ip,^Xgv[\'^n\'^ + |^n-iP]^n-i/2)- Using again (158|) . we can rewrite as in 
equation ( l59l) . where the terms cjn^ remain unchanged for all j = 1, 5, whereas the term Un^ 
now has the form 



(6) 



n 



, ,(6) , , /6) (6) 



where we use the same notation as introduced after (13H11 to define 

^n} ■= A5(y[|?/'„_i/2n(V'n-l/2 " U^^n + ^n-l)), 



For oOnl, expanding ipn-i and ■?/'„ around t = t„ — r/2 up to first order in r/2, we get , using 
(fTTj) and the mass conservation 



<ih2 < 4 l|A|Loo||\/|Loo(K.)-MN / dt Wmh^- (84) 



In order to estimate 00^21 we expand ipn-i and ■i/'n around t = t„ — r/2 up to first order in r/2 
and get similarly 



< a^'\a^:^T + a^:}^T' + a^^^r'), (85) 

where we define 

:= l|A|Loc||l^|Loc™.)-M[^o]'/', 



a' 



a« := M[^o] I dt 



L2 

^77 



^n^ ■= 2 ll^"-V2|lL2(||^n-l/2||i2 + / dt ||'(^^tj||^2 



<n-l 



al=^) := - Tdt ■■■■ ' 



t,i-i 
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For ujnl, using in particular again — \w\'^\ < \z + w\\z — w\ for all z,w G C and the 
decomposition ( l53l) . we get 

< ^l|Af7y[|^„|' + IVn-ll'll^n-^n + V'n-l -^n-l)|li2 

< \\M\L^\\y\\L^(R<i)iM^o]+M[llJoh]){\\pn\\ + \\pn-l\\ + PnW + ||^n-l||). (86) 

Therefore, plugging the estimates ( !84l) . ( !85l) . and ( !86|) into ( !60l) . we again find the closed inequal- 
ity the coefficients A^i^^ and A^^^ having the same form as in (ITOj) and (17T|1 . respectively. 
Using estimate flH^ on cj^^l, we see that the coefficient An'^ in the incoherent case contains all 
the terms from fl72|) of the coherent case with ci replaced by a^^-* , plus an additional term of the 
form a'^°)(an "* +an V + alfV^) which is due to the estimate of cj^^ Plugging the coefficients 

An\ An \ and A^^'' into the iterated bound fl7^ and using estimate (ITHjl on 6^0, we again get 
an estimate of the form ( 1771) where the constant C3 has the same form as in (1781) whereas the 
constant C4, compared to ( 1791) . now looks like 

C4 := |((a(°) + ||t;|L.)^'^dt ||#)|L. + ^^'^dt if^'WIL^ + ^'^dt ||A^(t)||^.) 

+C2a^'\M[M + |T max At 

max + / dt (87) 

Herewith, as in the proof of Theorem [31], we arrive at the assertion. □ 

Remark 38 Using estimates as in Remark [33], we can again bound the constants in (1H7|) explic- 
itly. 
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